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We consider quark-gluon plasma with chemical potential and study renormaliza- 
tion group flows of transport coefficients in the framework of gauge/gravity duality. 
We first study them using the flow equations and compare the results with hydrody- 
namic results by calculating the Green functions on the arbitrary slice. Two results 
match exactly. Transport coefficients at arbitrary scale is ontained by calculating hy- 
drodynamics Green functions. When either momentum or charge vanishes, transport 
coefficients decouple from each other. 

I. INTRODUCTION 

In the application of AdS/CFT [Tj, calculations are usually done on the holographic screen 
at infinity. However, according to the renormalization group (RG) ideas [2HI], those are the 
UV fixed point values which can not be reached by experiments performed at the finite 
energy scale. Therefore we need to run them down to the scale where one actually performs 
the experiment by the renormalization group flow. Recent studies of Wilsonian approach P, 
[7j to holographic RG flow of transport coefficients j5] in the framework of gauge/gravity 
duality suggest that those apparently different approaches of sliding membrane [9] , Wilsonian 
fluid/gravity [8] and holographic Wilsonian renormalization group (HWRG)P, [7J [TT] are 
equivalent [5]. Some of the transport coefficients such as shear viscosity 77, AC conductivity 
and diffusion constant D have non-trivial radial flows so that they interpolate the horizon 
values [UJ and the boundary values [H] smoothly. The holographic Wilsonian RG is also 
useful to understand the low energy effective theory for holographic liquid [121 H31 HE]- 

Since the discussions so far has been mostly for zero charge cases, one natural question 
is how to extend it to the case with finite chemical potentials. Unlike the zero charge case, 
metric fluctuation and Maxwell fluctuation will mix in charged black hole background. In 
general, the mixing effect of metric and Maxwell fluctuations in charged black hole is of 
essential importance, since it is the reason why transverse vector modes of Maxwell fields 
can diffuse and longitudinal Maxwell modes can have sound modes in the presence of the 
charge. Part of the answer has been given in [8] where the cutoff dependence of diffusion 
constant for the shear part of metric perturbations was calculated. The authors achieved it 
without explicit decoupling procedure by taking a specific scaling which, somehow, effectively 
decouples the mixed modes. However, it is not clear how to understand why it happens and 
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what scaling limit should be taken for other (i.g. sound) modes. 

In this paper, we first establish the flow equations of transport coefficients in the presence 
of chemical potential and numerically integrate the flow equations. We then calculate trans- 
port coefficients directly by hydrodynamic calculation on the holographic screen at finite 
radial position and compare the two results. We will find complete agreement. 

This paper is organized as follows. In section [XTJ we will briefly review the holographic RG 
flow of transport coefficients. We also derive the running diffusion constant in zero charged 
black hole from sliding membrane which reproduces the previous result in [SJ obtained by 
scaling method. In section III, we set up the charged black hole medium and focus on mixing 
structure of metric and Maxwell perturbations. We emphasize on organizing the coupled 
equations of motion. In section |IV[ we first write down the decoupled flow equation for 
both electric conductivity and "conductivity" for momentum current at zero momentum. 
And then we worked out the "conductivity" for momentum current in diffusion scaling limit 
following [8j. We write down the complete mixed flow equations for electric conductivity, 
momentum current "conductivity" and mixing parameter in the end of this section. In 
section [Vj we use the hydrodynamical method to calculate the Green functions at finite 
screen r = r c . We obtained the transport coefficients from Green functions by Kubo formula 
and we found complete agreement with the results from flow equations in section IV We 
conclude in last section. 



II. HOLOGRAPHIC RENORMALIZATION GROUP AND RUNNING 
TRANSPORT COEFFICIENTS: A REVIEW 

In this section, we shall discuss a few approaches to RG flow of transport coefficients and 
discuss the equivalence between them: They are sliding membrane paradigm [H] , Wilsonian 
fluid/gravity [5] and Holographic Wilsonian RG [EJ E]. The radial flow for transport coeffi- 
cients can be derived from holographic Wilsonian RG equation. It can also be derived from 
the classical equation of motion. The idea of membrane paradigm is to consider a mem- 
brane action pU] coming from the boundary action at sketched horizon. Similarly sliding 
membrane [9] paradigm uses the boundary action at arbitrary slice. From that action, one 
can obtain the retarded Green function by solving the equation of motion perturbatively in 
hydrodynamic regime. The linear response theory then gives the transport coefficients. The 
Green functions and transport coefficients satisfy the same radial evolution coming from 
holographic Wilsonian effective action Sb, which in turn can be obtained from integrating 
out UV geometry directly [5j [7] . Sb can also be treated as a boundary action at cutoff slice 
and it induces multi-trace deformations for IR dynamics [5J [7J . For more discussions about 
the essential multi-trace deformations see [231 121] • If one define the linear level transport 
coefficient at cut off slice, then the Hamilton- Jacobi equation for the effective action can 
directly give radial flow of them. These flows are exactly the same as those coming from 
equations of motion as explicitly shown in [5]. 
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Technically, sliding membrane paradigm is more convenient than holographic Wilsonian 
RG simply because the former starts from equation of motion while the latter is conceptually 
more satisfactory. We shall mostly use the former in this paper. 

We briefly review the sliding membrane paradigm [9J below. We derive the cutoff 
dependent diffusion constant, which is also calculated from the Wilsonian approach of 
Fluid/Gravity [§]. We show that the result from the sliding membrane paradigm repro- 
duces that in [8] to demonstrate the equivalence of sliding membrane paradigm [9] and 
Wilsonian approach of Fluid/Gravity [8]. We start with standard Maxwell action 
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with the background metric 



ds 2 = —g tt dt 2 + g rr dr 2 + g ii dx i dx l 



(1) 



(2) 



Defining J M and G by 

equations of motion can be written as 1 



G 



-gig 



eft > 



d t J l + d z J z = 
dr-J 1 + Gg u g zz d z F zt = 



d r J z - Gg tt g zz d t F : 



zt 







(3) 



(4) 
(5) 
(6) 



where z is the momentum direction and we focus on the longitudinal mode first. The Bianchi 
identities are given by: 



9rr9z 



-d t J z 



g rr gtt 



d z f + d r F zt = . 



G 1 G 

Using the definition of the conductivity a = J z /E z = J z /F zt , we have 

d r J z F zt - J z d r (F zt ) 2 (J z ) 2 



F 2 



a' 



(F zt ) 2 • 



(7) 



(8) 



We want to replace all the d r terms in above equation using equations of motion. Taking 
use of Q, ^ and (J7|, all the currents and fields disappear and finally we have the flow for 
conductivity [9] 



2 2 

9cs° 



-gg rr g zz 
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g 2 e s 



-gg tt g zz 



(9) 



1 Here we follow the notation in [9]. 
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Before doing explicit analysis for this conductivity flow we should mention that there are two 
important scaling limits for the frequency dependent flow equations: one is u ~ k 2 « 1 
and the other is A; = 0,<v << 1. In the diffusion scaling regime u ~ k 2 << 1, we have 
following solution of ^ for conductivity : 

i i i r- A 

rH 



a(r c ) oh iw/o f 
Since the (complex) conductivity at arbitrary slice and the green function is related by 

G R (kfj,, r c ) 



(J lJ (k^r c ) 
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we can write the Green function at cutoff surface in terms of horizon conductivity oh 



G l ^(k^, r c 



(jj 2 a H 



ioj — D(r c )k 2 



(11) 



(12) 



where D(r c ) = oh/ fo is the cutoff dependent diffusion constant. Substituting into (10), we 
have 2 

± = (15) 
a{r c ) o H ioj a H 

In the orthonormal frame, we have the normalized momentum as follows 

k 



CO — > 0J r 



k -> k c 



(16) 



Using these new variables we can rewrite (10) as 

1 1 k 2 D(r r 



o\r r 



where D(r c ) is given by 



D(r r 



D(r r 



9u 



(17) 



(18) 



Defining the local temperature by T c = , and a dimensionless diffusion constant by 



D c = DT C , we have 



Dc 



OH9U 

fo\/9tt 



(19) 



2 Remember that in this formula, the conductivity is defined using 



cr(z c ) 



J' 



-F t , 



For a physical observer hovering at z c surface we want to use the following normalizations 



a(z c ) &{z c ) 



(13) 



(14) 
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This is precisely the result obtained in [8], supporting the equivalence of two approaches. 
So far we considered only chargeless case. It is known that presence of charge introduces 
non-trivial mixing between modes. We will consider what happens for the flow of charged 
black hole below. 



III. MODE MIXING IN CHARGED ADS BLACK HOLE 



In order to describe low energy physics of various strongly coupled systems, we need dif- 
ferent IR deformations for the d+ 1 dimensional AdS background. Consider a d dimensional 
holographic system with finite charge density. The minimal d + 1 bulk action is 



S 



d d+l x 



2k 1 



R 



d{d-l) 
L 2 



(20) 



where Newton constant Gn = k 2 /87t and g is the Maxwell coupling. Maxwell and Einstein 
equations are given by 



-gFn 



. > \ -^V 2 



1 d(d-l) 



The charged black hole solution for this action is 
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ds 2 

L 2 K 2 



^.2 t 2 /-/'T'2 

I ,(-f(r)dt 2 + dx^) + - W) 



fJLil- 



„d-2 



f(r) = l + aQ 



2 r 



2d- 2 



^2d-2 



-(1 + qjQ 



2\ r 



(21) 
(22) 

(23) 
(24) 



Q 



(d-l)(d-2)g 2 



is a dimensionless coupling. Chemical potential /x is related to Q by 



and the horizon tq is the largest root of f(r) = 0. 



A. Modes in RN-AdS 



We consider the perturbations 



(25) 



around the background and A^, which are given in (23) and (24). One can use back- 
ground metric raise and lower tensor indices. The linearized gravity fluctuations can be 
decomposed into tensor, vector and scalar type jTHl US] . We will consider the first two types 
and leave the last one elsewhere. In the charged black hole background, vector part of metric 
perturbation and transverse Maxwell perturbation will couple. We can organize the the for- 
mer such that its action and equations look like that of longitudinal Maxwell perturbation. 
For more details about the action of linearized fluctuations in RN-AdS, see [2H 122] . 
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1. Tensor Mode bg x 

By symmetry argument or by direct analysis one can find that off-diagonal perturbation 
(f> = 5g x decouples from all other perturbations. Thus it satisfies the same equation of 
motion to that for a minimally coupled massless scalar: 

d,(V^gdy) = . (26) 

This evolution for will give the flow of correlation of the corresponding operator T%. Notice 
that if we work in d = 3, we should take k — > limit since we only have two spatial directions. 



2. Vector Modes b~g x z , Sgf , 5g x , 8A X 

Picking up a momentum parallel to z direction, we can only have d — 2 transverse di- 
mension left. One can observe that the shear mode of metric perturbations of components 
9xzi 9xti 9xr which behave similarly as Maxwell field a z ,a t ,a r [T7j defined by 

a z = 8g x z , a t ee 8g%, a r = 5g x r . (27) 

In terms of these variables, finally the vector part off shell action in d + 1 dimensional 
RN-AdS can be written as follows: 3 

8 = J d^x [~^F, V - Jg^rfr + j 2 a t A'^) ■ ( 28 ) 



Here we have used j and / to denote the current and strength for the a 



f = — YT^^-9f rtl > fi* = d » a v ~ d v a n , ( 29 ) 

9eS\ r ) 

and J and F to denote the current and strength for A^. Consider first the chargeless case, 
equations of motion for shear part metric fluctuations can be written as 

dtj 1 + d z f = (30) 

d r f + Gg a g zz d z f zt = (31) 

d r f - Gg a g zz d t f zt = , (32) 

where the effective determinant and the effective coupling are 

G = V^/9 2 cS (r), ^)=f5- ( 33 ) 



3 We derive this action up to total derivative terms, which are irrelevant to equations of motion. See 
Appendix A for details. 
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For the charged-AdS black hole, we introduce the new charge density which is defined as 

f = f - ^V^A' t A x . (34) 

The equations of motion takes the same form as chargeless case, 

dtf + d z f = (35) 
drj 1 + Gg tt g zz d z f zt = (36) 
d r f - Gg tt g zz d t f zt = , (37) 

but j l is replaced by j*, while the Bianchi identity holds as before 

- 9 ~^d t j z - 9 -^d z f + d r f zt = . (38) 

Defining J x and J x by 

J x = -- 2 ^~gg rr g xx d r A x , J x = r + - 2 V=gA' t a t (39) 
T T 

the equation of motion for A x can be written as 

- d r J x + j 2 ^~gg xx {-g u d t F tx + g zz d z F zx ) = . (40) 

3. Vector Mode A r at k 7 = 



From (40) and (36), one can see that A x decouples from a t in the k z — > limit: 

1 
9 



d r (^gg rr g xx d r A x ) - -^r gg l^ t ) 2 A x + u 2 ^-gg xx g tt A x = . (41) 



The relevant on-shell action for A x at boundary r = r c can be written as 

'S'on-shell = / J* A x . (42) 



IV. FLOWS OF TRANSPORT COEFFICIENTS 

We shall establish the transport coefficient flows for the charged black hole in this section. 

A. Shear Viscosity Flow and Scalar Response 

Since tensor mode decouples from all other perturbations, it behaves as a massless scalar 
perturbation even in the charged black hole background. Following the sliding membrane 
argument 0, it is useful to define a cutoff dependent scalar response function 

G,(r c ,K)= . (43) 
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And we define 

V(r c ,k,) (44) 
which satisfies a flow equation evolving equation of motion at k = 0: 

'2K 2 ri 2 (r c ,u) _ ^g_g u ' 

It is manifest from above equation that in the zero frequency limit, rj (r c ) is independent of 
r c . Its value is request to be 

V(r c ) = v{r ) = Y^{j) • (46) 



d rc 7]{r c ,uj) = iuj ( — — ) . (45) 



by the horizon regularity Since the entropy density is s = ^ (^) , the ratio rj/s does 
not run. In the next section we will see that this is consistent with the direct calculation of 
i](r c ) using the holographic hydrodynamics at the finite holographic screen. 4 



B. Electric conductivity flow at the zero momentum 



Now we shall consider the conductivity flow coming from A x perturbation with zero 
momentum. From (42) one can obtain the conjugate momentum for A x (r) at r = r c 



IL 



SS, 



on-shell 

5 A v . 



J x + 



-gA' t a t 



(47) 



Notice that this is nothing but the shifted current J x which was introduced previously 
to simplify the equation of motion. As explained in [19], the first and second term in right 



hand side of (47) should be related to electric conductivity and thermal-electric conductivity 
respectively. Notice that the on-shell action is the integral over the membrane at r = r c and 



the conjugate momentum eq.(47) is also defined at r c > r . In the limit of zero momentum, 



A x decouples and we can define electric conductivity by 

J x 

a A (u},r c ) 



iu>A 1 



One can quickly rewrite the equation of motion (|41|) as flow equation for 

d rc cr A 



+ 



-lUJ 



-ggrrgxx 



+ 



g 2 u 2 



9' 



-gg xx g tt 



. 



Due to the regularity condition at the horizon, we need 



cth = cTA(r ) 



1 



r \ d - 3 
T. 



(48) 



(49) 



(50) 



4 In reference [8], the entropy density is denned to be proportional to one over embedded volume of the 
cutoff membrane and r\ is also defined that rj/s does not run. As a consequence both their entropy density 
and viscosity run. Also both vanish at the infinity. 
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which is consistent with horizon conductivity evaluated in [9]. In order to see the solutions 
explicitly, we plot the conductivity flows with different black hole charges in d = 3 in Figure 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

1/r 1/r 



FIG. 1: r flow of AC conductivity, with d = 3 RN-AdS black hole background. Here we normalize 
r such that the horizon localizes at r = 1. 

A few remarks are in order for these flow solutions. 

• Zero Charge Limit: Remember that the charge density in the boundary theory is 
related to the chemical potential by 

p = pr d - 2 /9 2 ■ (51) 

When /i = chargeless case, the d = 3 conductivity is independent of the cutoff 
position. The flow solution for d = 3 is trivial while there are nontrivial flow solutions 
for d ^ 3 [5]. When d ^ 3 there is a fixed point at boundary but not at horizon. 
This is because near the boundary metric is aysmptotically AdS and scale invariance 
is there. This scale invariance is lost when two fluctuating modes mix in the RN-AdS 
black hole case. 

• Extremal Limit: As shown in Figure [TJ when the charged AdS black hole becomes 
extremal, one can find a fixed point from flow of conductivity near the horizon due 
to the appearance of AdS2 near the horizon. This fixed point will disappear in non- 
extremal case. One should notice that there is no fixed point for the conductivity near 
the boundary both for extremal and non-extremal case. The evolution equation for 
a a loses scale invariance near the infinite boundary in the presence of charge due to 
the mix of at, A x modes. 

• Check against Boundary Result: In order to check the consistence with previous results 
calculated at infinite boundary r = oo, we plot the boundary AC conductivity in Figure 
[2j which is precisely consistent with the results in |l9j . The Q dependent term in the 
flow equation clearly give the origin for the divergence of the imaginary part of the 
conductivity in the zero frequency limit. 
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• Conductivity Minimum: It is interesting to observe that there is a window of pa- 
rameters u and fi where the real part of conductivity curve in radial direction has a 
minimum as shown in Figure [TJ The minimum of the conductivity pick up a certain 
scale r* and we show the charge dependence of that scale in Figure [2j 




1 2 3 4 5 0.0 0.5 1.0 1.5 

ai f/ A 2 



FIG. 2: Left: Frequency dependence of boundary AC conductivity, with d = 3 RN-AdS black 
hole background . As a check for our flow solutions, the boundary AC conductivity is precisely 
consistent with the results in [19]. Right: Chemical potential dependence of scale r* with minimal 
conductivity in d = 3 RN-AdS, where \x is dimensionless chemical potential. 



C. RG flow of 'Conductivity' of momentum current flow at zero momentum 

Since the metric perturbation a M mimic the Maxwell system, we can define the "conduc- 
tivity" for the mode a z = h z by 

<?h ■■= T~ ■ (52) 

Jzt 

At k z = 0, since a M modes are decoupled with A x , we can obtain the decoupled flow of ah 
as follow 

d r O~h 2 9rr9zz n J 1 r . zz n /cr 

: + Oh G 9 9 = • (53) 

This flow equation has the same form as the longitudinal conductivity flow in zero charge 
case although the metric components contain Q dependence. As a consistent check, this flow 



equation is exactly the same as shear viscosity flow in (45 ), because when k z = 0, there is no 
polarization direction and ah is nothing but shear viscosity by definition. We will explain 
more about the physical mean of ah in next subsection. 



D. RG flow of 'Conductivity' of momentum current in diffusion region 



Since shear parts of metric perturbations behaves as longitudinal Maxwell field with a r 
dependent effective coupling, we shall work out the "conductivity" flow in diffusion region. 
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In order to handle the equations of motion: (35)~(38), we write them in momentum space 
and take the diffusion scaling following j8] 



t 



(54) 
(55) 



First consider the in-falling boundary condition, which requires j z linearly related to f zt 
near the horizon. This condition requires us to scale j z as same power of f zt : 



r 



e 3 f(0) +ef (D + 



Through the charge conservation equation we obtain 



and we can also obtain 



e [r +ej + 



(56) 



(57) 



(58) 



since there is no different scaling between j f and j t . Now we want to take e — > limit, and 
only the lowest orders of fields are left. From (36) and (37) we obtain 



3 r f m = 



8,f" s > = . 



Requiring charge conservation, we have 

f 



(o)=J t (°)^ 



From (38), we have 



9r/ (0) = 



G 



It means we have the solution for / 



(0) 



zi 



f2\r) = f£>(r )+2K 



(0), 



g 2 d z j t(0 \r) 
s/^99 rr 9 a 9xx 



dr 



(59) 

(60) 
(61) 

(62) 



Following the usual definitions in sliding membrane for chargeless case, we define the con- 
ductivity by current and electric field as 



J zW (r c ) _ k f (0 \r c ) 
f$\r e ) <of£\r e ) 



Remember 



r t (o) 

J =3 



(63) 



(64) 



is a constant. By taking the scaling limit we have the simplified equation of motion for A 

1 



(o). 



g 2 



(65) 
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The solution for can be solved explicitly as shown in [8]. 

Since j z ^ is a constant in radial direction, after having the solution for (|62|) can be 
rewritten as 

k 2 D h (r c )/T c 



1 

crh(r c 



1 

c-fc(ro) 



iuj c a h (r ) 



where the horizon value is 



1 



"[0 

16nG N \L 



d-i 



(66) 



(67) 



This is exactly the "conductivity" flow for shear part of metric fluctuations. This is our main 
result in this subsection. After taking the scaling limit following [8J, we have the analytical 
result for RG flow of "conductivity" of momentum current. 
As a byproduct, we obtain the D h as 



D h {r c ) 



1 



(rf-2)(l-«Q 2 ) + 2 



4tt (1 + aQ 2 )d - 2(d - l)aQ 2 (r /r c ) d - 2 



(68) 



which is charge dependent diffusion constant. Notice that we have used the solution for A x 
to obtain the Dh- This result was first obtained in [8], where the diffusion constant is derived 
from Fick's law, while here we find that this diffusion constant is included in conductivity 



flow for momentum current without any boundary condition for in (|62|). 



Now we will explain the physical mean of the "conductivity" ov Notice that a z = 5g x 
corresponds to boundary operator T*. At zero frequency retarded Green function of T* is 
exactly the same as T% since there is no special polarization direction. Thus the transport 
coefficient given by a z = 5g x will correspond to longitudinal momentum dependent viscosity. 
This can be confirmed at the horizon, since o^, at horizon is nothing but [9] 



s 

Air 



which equals to shear viscosity. This is consistent with rj/s 
transport coefficients at the horizon are all frequency independent 



j_ 



(69) 



One should note that 



E. Mixed RG Flow Equations 

We shall discuss the mixed flows coming from coupled equations of motion without any 



scaling limit. Consider again the coupled equations of motion: (35) - (38) and (40). By 
defining 

3 ^ ^ 
Oh = -r- , °a = —-r- , a := — (70) 
fzt tuAz A x 



we can derive the following flow equation for them as follows 



d r a h 
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Gg tt g 



(71) 
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d r a A 



—iu 
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-9K 



-99 



-9 U + 9 Z 
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—iuj 



tt „zz 



Gg tt g 

Oh 



9 2 °a 



-gg xx g rr 



o, (72) 
(73) 



One can easily observe that when Q = 0, all mixing effects disappear and (71) and (72) 
will reduce to longitudinal and transverse form of conductivity flow for chargeless case [S] 
respectively. Eqs. (71), (72) and (73) give the mixed RG flow. 



V. HYDRODYNAMICS AT THE FINITE HOLOGRAPHIC SCREEN 

In this section, we study the RG flow of transport coefficients by calculating the Green 
functions at finite screen r = r c . We consider the shear mode in the 5 dimensional RN-AdS 
background. This hydrodynamics problem is considered in [21 J for r c — > oo case. We use 
the following gauge 

a r = , A r = (74) 
and introduce following coordinate and notations: 

u = ^ , a = aQ\ b = — . (75) 
r z 2r 

In terms of the rescaled vector B(u) = A x (u)/fi, the equations of motion are given by 

lb 2 / \ 

= a t " a t ' [ojka z + k 2 a t ) — 3auB' , (76a) 

u uf \ J 

= kfaj + ua t ' — 3auuB, (76b) 
= a z " + K —^aJ + — (u 2 a z + uka t , (76c) 

= B" + ^B'+^-L 2 - k 2 f) B - lot'. (76d) 

Notice that B and a t couples which is the source of the complication. To handle the problem, 
we introduce the master fields [25J, 

$ ± = -a/ - 3aB + —B, (77) 
u u 

with C± given by 



C± = (1 + a) ± 7(1 + a) 2 + 3ab 2 k 2 . 
The decoupled differential equations satisfied by the master fields are 

(« 2 /)' , , b 2 /, ta x C, 



= $ ± " + ^U±' + ^ - fc2 /) $ ± - (78) 
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We expand the master fields in the hydrodynamic regime in following way: E 



$+ = g(u)C(l + buoF +l {u) + b 2 k 2 G +2 {u) + b 2 uo 2 F +2 (u) + 
$_ = (7(1 + buF^{u) + b 2 k 2 G^ 2 {u)b 2 u 2 F^ 2 {u) + •••), 



(79) 
(80) 



where the factor singular near the boundary is explicitly taken out for $ + and it is given by 



9W 



3a 



u 2(1 + a) 



(81) 



In this regime, the equations of motion are solved by imposing the ingoing boundary condi- 
tion at the horizon in [2T]. Here, we collect a few terms which will be relevant later 



KM 

FU(u) 



G'_ 2 {u) 



(2 -a) 



4(1 + a) 2 u 2 f(u)g 2 (u)' 



u 2 f(u) 



2(1 + a)u 



2 ' 



(82) 
(83) 
(84) 



The constants C and C are integration constants and fixed by imposing the boundary 
conditions. It is convenient to define a gauge invariant field Z by 



ka t (u) + ua z {u) = Z. 
The gauge field B{u) can be expressed in terms of master fields as 



B 



1 



Using (76a) and (77), we obtain 



u 2 <$>' ± - C±uB' = jkZ - C±B. 



(85) 



(86) 



(87) 



With (86), l.h.s. of (87) can be expressed in terms of master fields Q±. Requiring (87) at 



u = u c we can determine C and C in terms of the boundary values of Z and B at u — u, 

a(u c )b 2 kZ c - /3(u c )f(u c )B c 



C 
C 



u c g(u c )(-iu + D{u c )k 2 )b 
a(u c )b 2 kZ c - (3{u c )f{u c )B c 
u c g(u c )(-iu + D(u c )k 2 )b 



(88) 
(89) 



Here, the definition of F +1 (u) is related to the Fi(u) in [5T] as 



bF +1 (u) =F 1 (u) 



4ttT 



log(l - u), 



and similarly for F_i(u), etc. 
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where denominators are expressions up to 0(u> 2 ) and 0(k 3 ), and D(u c ) is given by 

D( u ) = = 

2(1 + a)u c g(u c ) 2(1 + a) — 3au c 

The coefficients a, (3, a and (5 are given in terms of the solutions at u = u c as 

a{u c ) = -u c g{u c ) (l + buF +1 (u c ) + b 2 k 2 G +2 {u c )) H 

„, , 3a6 2 fc 2 

/3(Mc) = 2(rT^) + --- 

fe 2 £; 2 

a(« c ) = -ibu + — + ■ ■ ■ 
2(1 + a) 

5(u c ) = 2(1 + a)ibu - b 2 k 2 



(90) 



(91) 
(92) 

(93) 
(94) 



Now we calculate the Green functions at u — u c . We start from the Einstein-Hilbert 
action with Gibbons-Hawking terms and counter terms [26] for the gravity part 6 



ravity 



Seh + Sg.h + S, 



c.t 



(95) 



where the Gibbons-Hawking terms and counter terms can be expressed in terms of (trace 
of) the extrinsic curvature K and induced metric 7 MJy as 



1 



16nG 



d d+1 x 



i 



gR 



c.t 



8ttG 



(96) 



Both Gibbons-Hawking term and counter terms are defined as the hypersurface at u = u c . 
After perturbing the action and integrating out the classic solution of perturbations, finally 
we obtain the boundary action for the shear modes, which are given by 



5° 



gravity 32^4 



d 4 x 



1 , 3 

-a t a t 

u u z 



1 



1 



a t a t 



M aza > + L U f{u) _ 3 Jj@j _ u f {u) \ 
u u z V / 

d 4 x [-3af{u)BB' + 3aBa t ] . 



a,a. 



(97) 
(98) 



-gauge 32^2^4 

Once we have C , C in terms of boundary values, by the definitions of master fields we can 



6 We use this action in order to keep consistence with previous results when u c — > 0. 
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express first derivatives of a t , a z and B in terms of boundary values at u = u c as follows: 

// x _ at(u c )b 2 kZ c - /3 t (u c )f(u c )B c 
at[Uc) ~ u c g(u c )(-iuj + D{u c )W)b ^ 

// \ _ a z (u c )b 2 kZ c - /3 z (u c )f(u c )B c 
aAUc) ~ u c g{u c ){-iu + D(u c )k*)b UUUj 
R // x a B (u c )b 2 kZ c - /3 B (u c )f(u c )B c 

[Uc) u c g{u c ){-iu + D(u e )k?)b ' 1 J 

where the coefficients are explicitly given by 

a t (u c ) = -u 2 c g{u c ) [I + bu(F +1 (u c ) + F_i(u c ))) + 6 2 fc 2 (G +2 (« c ) + GL 2 (u c ))] + • • • , (102) 
^ / \ 3aibui 

Pb{u c ) = - + • • • , 103 

2(l + a)/(« c ) 

where . . . denotes higher frequency and high momentum terms. Notice that we only write 
down two coefficients in above equations for later use and leave other coefficients in Appendix 
B. With the help of above results, the Green functions at slice u = u c can be read off from 
on shell action 

L 3 / at(u c )b 2 k 2 3 / 1 \\ 

Uxtxt 16 K W\ulg(u c )(-iuj + D(u c )kZ)b ul{ yfj^)]' [ } 

G xzxz = L ( Oi t {u c )b uj + J_ faf( u j _ s^f( Uc ) - u c f'(u c j) X , 

16/sro 4 yu z c g\u c ) (— iui + D(u c )k z ) b u z c V /J 

(105) 

r _Jl 2a t {u c )b 2 ku 

xtxz 16/€ 2 6 4 ^( Mc )(_^ + jD ( Mc )A; 2 )6' 

r L 3 \ -^{u c )f{u c )-3au c f{u c )a B {u c )b 2 k 2 \ 

xtx 32« 2 &Vl ~ ^ 5 K)(-ia; + J DK)A; 2 )& J ' 1 J 

r £ 3 /3 z (u c )f 2 (u c ) -3au c f(u c )a B b 2 ku 

xzx 32k 2 6V «^(« c )(-iw + D(u c )k 2 )b ' 1 j 

G _ £ 3 3a/ 2 K)^K) 
x * 16k 2 6V^K)(-^ + ^K)A; 2 )6 ' 1 ; 

A. Cut off dependence of diffusion constant 

One can easily observe a universal diffusion constant depending on the cutoff position 
from all the Green functions: 

D(u c ) = - bfM , , = - hf ^\ . (110) 
v 1 2(1 + a)u c g(u c ) 2(1 + a) - 3au c v 1 

Change to orthonormal frame, one can obtain 

D(u c ) = D(u c )^L . (Ill) 
V9tt 
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Together with a normalization with local temperature, one can obtain the dimensionless 
diffusion constant 



D c = D(u c )T c = D(u c )T H 
1 2 - a 



9zz 

9tt 



4vr2(l + a) - 3au 



(112) 



This is nothing but in (68). The above result was first obtained in [8] by taking a certain 
scaling for the equations of motion. Here we show that this diffusion pole appears in all the 
cut-off dependent Green functions. 

B. Shear viscosity 

The shear viscosity is calculated by using Kubo formula 

. . , lmG xvxv (cj, k — 0, r c ) 
77(77.) = - lim y yy - 113 

For k = 0, z-direction can be treated equivalently to y-direction, since there is no polarization 
direction, and we have 

L 3 

G xyxy (u, k = 0, r c ) = G xzxz (u, k = 0, r c ) = ~ i ^^ UJ + 0(u 2 ). (114) 



Then the shear viscosity becomes 
which is constant independent of cut-off. This is consistent with the result from flow equation 



(46). 



C. Cut off dependence of DC conductivity 

The conductivity can be calculated form the O(oj) term of the Green function G xx . 
Generally it contains the 0(u 2 ) terms of the master fields and have a complicated expression. 
Taking into account the 0(u) 2 ) terms, it becomes 

q _ _Jl_ 3af 2 (u c )P B (u c ) 

xx IGK^^UcgiuJ^itu + Diu^ki-ulfiuJiF^ + F^FL!))^ 1 ' 

where /3b{u c ) also have the 0(u> 2 ) corrections as 

p B (u c ) = ~ 2(1 ^ K) + bW [(ug(u))' FL 2 (u) + (ug(u)F +1 )' FL^u)] (117) 
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For k — > 0, the expression of G xx is simplified and given by 



G xx = °-^££r ( ^Ff + ^F' +X {u c ) ) + C^ 2 



3af(u c )L 3 ( (u c g(u c ) )' 
lQK 2 b 4 n 2 V u c g(u c ) 

3af(u c )L* f M(« c ))' + (2 -«) a 1 . + 0(a/l) (ng) 



16/t 2 6V V ^K) 4(l + a) 2 M c 2 /K)(7 2 K) 

The real part electric conductivity is explicitly given by Kubo formula 



lr (2-q) 2 1 



where (7 is the gauge coupling. At the horizon, a DC is 

1 r 1 7rL 



°h = - 2T = - 2 ^(T+ v/T 2 + /i 2 /37r 2 ), (120) 
5- 2 2 

which is consistent with eq.(58) in [U] in the limit where charge or chemical potential goes 
to 0. This horizon value is related to the membrane conductivity a H = .P mb /E l [9] by 

V 9tt9rr 1 M s 

ctr = — j^^gu- a H = —. (121) 

One can also check that as one goes to the boundary {u c — > 0), DC conductivity is reduced 
to 

1 ro (2 - a) 2 
g 2 L 4(1 + a) 2 



which precisely agrees with previous result in 

Finally we can check the consistency of flow equation of AC conductivity by comparing 
its numerical value in the limit of zero frequency with that of DC conductivity calculated 
here. The result is plotted in figure [3j 



VI. CONCLUSION AND DISCUSSION 

In this paper, we present a method to work out the RG flows for transport coefficients 
for quark gluon plasma at finite chemical potential with charged AdS black hole dual. Due 
to the mixing effect between Maxwell and metric perturbations, we need to solve the cou- 
pled equations of motion, which is usually difficult. We organize the system as two coupled 
Maxwell systems and define two conductivities for each of them. With a parameter charac- 
terizing the mixing effect, we write down the mixed flow equations. These mixed RG flow 
equations will be simplified in certain limits. These flow equations will characterize how the 
transport coefficients will change as energy scale changes. 

We explicitly give the flow equations for conductivity and shear viscosity. In order to 
check these results analytically we use hydrodynamic method to calculate the Green function 
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Dashed: AC flow at infinite small frequency 
Solid: DC flow 
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FIG. 3: Checking consistency between flows of DC conductivities: Left: Plot of equation (119) 



from hydrodynamics. Right: comparing the hydrodynamic result and flow equation result. We 
can see they are explicitly same with each other. 



at finite cut-off slice r = r c . We impose equations of motion at r 



which is guaranteed 



by RG invariance of bulk action at classical level [5]. Then the Green function can be read 
off from the on-shell action at r = r c . We extended the usual counter term to arbitrary 
slice in order to have consistent result when r c — > oo. By using Kubo formula we obtain 
the analytical results of RG flow formulas for transport coefficients and we found complete 
agreements with that from flow equations. 
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Appendix A: Derivation of bulk action 

Starting from Einstein-Hiibert action, we calculate the bulk action for vector perturba- 
tions a z ,a t ,a r , A x in the 5D RN-AdS background. 



S 



d 5 x 



A l u' l (/ / + A 2 n[(/ -rDA' r A' r + E<i l A' r + [ -B iat a t 



1 



-B 2 a z a z + 



• (Al) 
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Here all the prime denotes the r derivative and ••• represent the terms with dt,d z . The 
coefficients are given by 

Ai = —z — :r, A 2 = „ = (A2) 
2k 2 2L 5 ' 2k 2 2L 5 v ; 

1 r 4 / rf'\ 1 2r 4 / /A . 

D - 1 " r3/ F 1 — 1 ^ fA4^ 

D " ?w (A4) 

Notice that the total derivative terms can be deleted for the purpose of the equations of 
motion. One can evaluate equations of motion from the action for a t , a z , A x : 

A x a" + A[a! t + EA' X = (A5) 
A 2 a" + A' 2 a' z = (A6) 
(DA' X + Eat)' = 0. (A7) 

However we should calculate the final on-shell action by adding proper local counter terms, 
as we did in section |V} We conclude that the mixed bulk action can be written equivalently 
as two Maxwell actions plus the only mixing term EatA' x . 



Appendix B: Coefficients for a' t , a' z and B' 
The coefficients for a' t , a' z and B 1 are given by 

a t {u c ) = -u 2 c g{u c ) [l + bu{F +1 {u c ) + F_i(w c ))) + b 2 k 2 {G +2 {u c ) + GL 2 (u c ))] + • • • , (Bl) 
<y z {u c ) = - rjr r a t (u c ) (B2) 

, . u 2 g\u c ) ^ (sW W - ^+1 W) - (f-i K) + i'+i W)</W) 

oc B {u c ) = — r + 



+ 



2(1 + a) 2(a + l) 

b 2 k 2 u 2 c (g(u c ) (G'_ 2 (u c ) - G' +2 (u c )) - (G- 2 (u c ) + G +2 (u c ))g'(u c )) 

2(a + l) 



3au 2 c g'{u c )b 2 k 2 

+ 4(l + a)3 + '--' (B3) 

/ x ., 3au 2 q(u c ) /rt . 

Pt(u c ) = ibco f\ c) +■■■ , B4 

f{Uc) 

n . . 3au c b 2 ku 

w " J = -^) + -' (B5) 

2(1 + a)/(tt c ) 
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